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1 Introduction. 



Quantum chromodynamics, (QCD), is widely accepted as the quantum field theory of the nucleon 
partons known as quarks and gluons. At high energy, perturbation theory provides an excellent 
description of the phenomenology of hadron physics from, say, the point of view of deep inelastic 
scattering. Although the physics of QCD at low energies is not as well understood, it is generally 
accepted that quarks are confined and the perturbative vacuum, which is used for higher energy 
QCD calculations, is unstable, ^|2J|3J|3]. To probe the infrared regime from the field theoretic 
point of view several formalisms have been developed. One in particular makes use of the 
operator product expansion and sum rules where the effect of dimension four operators such as 
(G^,) 2 are incorporated, jSJ- Such operators have non-zero vacuum expectation values and can 
be regarded as probing the structure of the true vacuum. Such methods have been successful 
in revealing insights into the low energy structure of Yang-Mills theory and QCD. However, 
in this general context there has been a large amount of activity recently into the effects that 
dimension two operators have on this scenario, |H1 d El El E01 HH [T^l [IS1 H31 HH1 EHl CZl [TH1 HH1 
For instance, the operator \A a2 and the related BRST invariant operator \A a2 + ac a c a have 
received particular attention being motivated by early work on the Curci- Ferrari model, |2(J| . 
where is the gluon field and c a and c a are the ghost and anti-ghost fields of the covariant 
gauge fixing which is parametrized by a. The BRST algebraic properties of these operators has 
been extensively studied, [HI E31 EH HS1 E3 H3 HE]- particularly in relation to unitarity. 
Unlike (G® u ) 2 such operators are clearly not gauge invariant and therefore they can only arise 
in gauge variant quantities such as the strong coupling constant. (See, for example, [29} 130].) 
However, one can construct a dimension two gauge invariant operator which is non-local, 
but is believed to reduce to \A a2 in the Landau gauge where it is clearly local. Therefore, 
various groups have examined the effect such an operator has in Yang-Mills theories both in 
the case of ordinary linear covariant gauge fixing and in a related non-linear covariant gauge 
fixing. The latter is known as the Curci-Ferrari gauge having its origin in the Curci-Ferrari 
model which was an early attempt to construct a field theory where the gluon and ghosts were 
classically massive, |2()| . This gauge differs from the usual covariant gauge fixing due to the 
presence of a renormalizable four-ghost interaction and a different ghost-gluon interaction. Its 
ultra-violet behaviour is equivalent to that of the usual linear Landau gauge when the covariant 
gauge parameter is nullified. 

One important study of the effects of the \A a2 operator has been provided in I31j . 
There the effective potential of this operator has been constructed at two loops in Yang-Mills 
theory. This is not a straightforward exercise as one has to apply the local composite operator, 
(LCO), formalism which correctly accounts for a composite operator being included in the path 
integral formalism, [111 1311 IH2l I33| . Early work with this technique centred on studying mass 
generation in the two dimensional Gross-Neveu model where estimates for the mass gap were 
obtained which were in good agreement with the known exact mass gap of the model, |32| I33j . 
Consequently the LCO technique provides an effective action which corresponds to the Yang- 
Mills action as well as a new contribution which involves an extra scalar field whose elimination 
by its equation of motion is related to the dimension two composite operator. In this respect it is 
similar to the Gross-Neveu model though in that case the analogous scalar field arises naturally 
in the Lagrangian. With the new action derived with the LCO technique the effective potential 
of the dimension two operator has been constructed at two loops in Yang-Mills theory, |llj . 
However, as the nucleon world possesses quarks there is clearly a need to repeat the analysis for 
full QCD. This is the main purpose of this article where we will extend the effective potential 
of to include massless quarks. Whilst it may appear that our calculation follows it 
should be noted that considered only the SU(N C ) colour groups and, moreover, required 
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the evaluation of three loop massive vacuum bubble graphs using the tensor correction method 
developed in By contrast, we will consider arbitrary colour groups which allows one to 

study, in principle, the effect in grand unified theories. More importantly, though, we will 
bypass the need to compute three loop massive vacuum bubble diagrams. These were required 
to determine the explicit expression for a particular parameter, C(sO> as a function of the coupling 
constant, g. This quantity, which plays a role similar to the coupling constant for the new scalar 
field of the new action, obeys a particular renormalization group equation and its explicit value, 
as we will demonstrate, only requires the evaluation of three loop massless Feynman diagrams. 
The importance of this observation is not to be underestimated since from a calculational point 
of view it reduces the work significantly and will be important in other similar applications of 
the LCO formalism for which the symbolic manipulation programmes developed here can be 
applied. One final motivation for this article is to ascertain the effect the quarks have on an 
effective gluon mass. In the construction of the two loop effective potential it transpires that 
for both Yang-Mills and QCD the naive perturbative vacuum is not stable, but there is a 
stable vacuum where the gluon field develops a mass through the non-zero vacuum expectation 
value of the new scalar field arising from the LCO formalism. Clearly an estimate of such a mass 
is important when quarks are present. While our calculations are still effectively perturbation 
theory probing towards the true non-perturbative vacuum, the full vacuum expectation value 
of the scalar field will be comprised of two parts. One is the perturbative piece estimated 
from the effective potential but the other component will arise from purely non-perturbative 
phenomena such as instantons, In this context there has been various studies of dimension 
two operators and condensates on the lattice giving estimates of similar quantities, |29l 130] . 
Hence, it will be useful to provide information with quarks present ahead of lattice simulations 
of full QCD where it would be hoped the dimension two condensates could be measured at 
an energy scale compatible with the energy scale we will determine our effective gluon mass 
estimates. 

The paper is organised as follows. In section two we will review the construction of the LCO 
formalism for full QCD and determine several quantities which are central to the construction 
of the two loop effective potential in the QCD case. The full two loop potential is discussed 
in section three whilst we provide a detailed analysis of the colour groups SU(2) and SU(3) in 
section four. Concluding remarks are given in section five. 

2 Formalism. 

In this section we review the basic formalism for constructing the effective potential of a com- 
posite operator using the LCO formalism, Our notation parallels that of and we will 
focus on the main ingredients as well as where new features emerge. The starting point is the 
construction of an effective action which includes the composite operator we are interested in, 
\A a2 , in the Landau gauge, and is based on the usual QCD action. In particular the QCD 
Lagrangian is 

L = — X -Gl v G a » v - -^-(d^A*) 2 - ^d^D^c a + ijPflp 1 (2.1) 

where Ct v = d^A^ — d v A^ — gf abc A h ^A1, f abc are the colour group structure constants, tjf 1 is 
the quark field and 1 < a < Na, 1 < I < Np and 1 < i < Nf with Np and Na the dimensions of 
the fundamental and adjoint representations respectively and Nf is the number of quark flavours. 
The covariant derivatives are defined by 

ZV a = ~ 9f ahc <c c , Drf" = d^ u + igT a AW J . (2.2) 
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Although most of the recent activity in this area has dwelt on Yang-Mills in the Curci-Ferrari 
gauge since it has properties in common with the maximal abelian gauge, we have chosen to fix 
the gauge with the usual linear covariant gauge fixing. This is because the potential in either this 
gauge or the Curci-Ferrari gauge is the same when a = 0. For reasons which will become apparent 
later the quarks in 1)2. 1|) are massless. Moreover, we can regard (|2.1|) as the bare Lagrangian 
and introduce renormalized quantities via the usual definitions such as Aq^ = \fZ~A A afl and 
g = Z g g where the subscript, D , denotes a bare quantity. 

To include the composite operator, \A a2 , in the path integral it would initially seem that 
one should include the term \ J A^ 2 in the Lagrangian where J is a source term and derive the 
usual path integral generating functional, W[J]. Ordinarily in such an approach one couples the 
fields themselves to each source to obtain a generating functional which is used to determine 
the Green's functions. However, in coupling to an object quadratic in the fields an immediate 
problem arises in that the new action ceases to be multiplicatively renormalizable, This 
is due to the generation of divergences proportional to J 2 which are related to vacuum energy 
divergences, jllj . Divergences involving one power of J are absorbed into the renormalization 
constant for the renormalization of the composite operator, \A a2 , itself, j31U34j . Similar features 
are also present in the Gross- Neveu model, (321 EH]- To circumvent this lack of multiplicative 
renormalizability one must add a term quadratic in J with an appropriate counterterm in order 
to have a sensible generating functional, W[J]. Therefore, the generating functional which is 
the starting point of the LCO formalism is, jllj . 



-W[J] 



VA^V^ Vil) Vc Vc exp 



d d x 



L, 



- T A a2 
2 J o /± on 



+ 



-C J 2 
2 Co J Q 



(2.3) 



where all quantities are bare and the new quantity Co has been introduced to ensure one will 
have a homogeneous renormalization group equation for W[J], whose conventions we follow 
throughout. 

With the extra terms in the action additional renormalization constants are required over 
and above those needed to render the normal QCD Lagrangian, L, finite, However, it is 

important to realise that the additional terms do not affect the explicit values of the ordinary 
QCD renormalization group functions. This can be understood from the fact that the new 
cubic interaction can be regarded as a mass term for the gluon which, from the point of view of 
renormalization, will regularize infrared infinities and leave the ultraviolet structure unchanged. 
Therefore, we need to introduce two new renormalization constants which to be consistent with 
jllj we denote by Z m and 5^ where the latter is strictly a counterterm. Therefore, using 



Jo 



Q J 2 Q = (c + so J 2 



(j2.3|) becomes 



-W[J] 



VA^V^V^VcVc exp 



Jd d x(l - \z m JAf + I( C + ,5C)J 2 ) 



(2.4) 



(2.5) 



For the construction of the effective potential the explicit values of each renormalization constant 
will be required. It turns out that these are straightforward to compute. For instance, Z m has 
already been determined in j2f)J ESI EH EH EHj at one and two loops and to three loops in |3I] . 
The method to achieve this is twofold. One can either use a Lagrangian with a gluon of squared 
mass J and renormalize with a massive gluon propagator in the Landau gauge or one can regard 
the operator ij A a2 as a composite operator and renormalize it as an operator insertion in a 
gluon two-point function using massless propagators. Indeed it was the latter approach which 
was used to deduce the three loop operator anomalous dimension [31]. However, a key property 
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that emerged from that work, which has subsequently been proved to all orders in is that 
the renormalization of is not independent. In the Landau gauge the anomalous dimension 
is given by the sum of the gluon and ghost wave function anomalous dimensions. As this will 
play an important role in the construction of the effective potential we note that in the Landau 
gauge the explicit value in the MS scheme is 



7m (5) I 



= 



35 
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+ 
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V 75607 



T F N f 
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C A T F N f 
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C\ + Y T F N f°F + -^ T F N f C A + ±T F N f C F 
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24C(3) + ^)7/.V / CVC' 1 



9 



(16vr 2 ) 3 
(2.6) 



ab 



where the Casimirs are defined by T a T a = C F I, f acd f bcd = C A 5 ab and Tr (V a T fe ) = T F 5' 

Therefore all that remains is to evaluate the counterterm 6£. In the source, J, was 
considered to be constant giving the gluon a mass. The counterterm was then deduced by 
renormalizing three loop massive vacuum bubble diagrams generated from the interactions of 
L. This required the development of the tensor correction method, [35 ; to efficiently handle the 
tensor reduction and algebra associated with massive Feynman diagrams. Then the divergences 
of the vacuum diagrams, which were proportional to J 2 , were absorbed into 6£. In our approach 
we have considered the source, J, to be a field which interacts with the gluon leaving no candidate 
for a mass term for the gluon. It therefore remains massless. By viewing the source like an 
interaction rather than a mass, we have avoided the need to evaluate massive vacuum bubbles. 
Instead we proceed by considering the vacuum Green's function 



(0\O(x)O(y)\0) 



(2.7) 



\J . It is the divergences of (j2,7|) which must be cancelled by 



using (|2.1|) where 0(x) 

5C,. We carried out the three loop renormalization by using the Mincer algorithm, jlU], as 
implemented in Form, |41l I42j . with the relevant Feynman diagrams generated by Qgraf, |43j . 
The converter files mapping the Qgraf output into input Form notation have already been 
used in [33|. However, the Mincer algorithm evaluates the divergences of two-point functions 
whereas our vacuum Green's function has no external legs. It was therefore necessary to cut 
each diagram through the one internal J propagator and consider the resulting J two-point 
function. This is valid because the propagator of this source field is a constant and is thus 
common to each graph. It is important to note that this is not equivalent to considering the 
two-point source Green's function for the complete Lagrangian, (|2.5|) . as this would contain extra 
diagrams with internal J propagators which do not contribute to (|2.7jl . From a calculational 
perspective it would seem that renormalizing (|2.7|) in this way is more efficient since one needs 
only to use massless propagators. Indeed this alternative could probably be applied to similar 
computations of the effective potentials of other dimension two composite operators. One benefit, 
for instance, of using the Mincer algorithm, |40| I41j. is that the tedious tensor reductions are 
automatically implemented. To check that this interpretation is consistent with the three loop 
massive vacuum bubble calculation of 13 lj we have renormalized (|2.7|) for arbitrary Nf and 
found in dimensional regularization with d = 4 — 2e using MS, that 
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+ \™T*Nf - 4T F N f C F - ^-T F N f C A + ^ 
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When Nf = (|2.8|) corresponds exactly with the result of [TT]*. 

With (|2.8|) we need to construct the associated renormalization group function denoted by 
5(g) and defined by 



dc 



dfj, 



2jm(g)( + 6(g) . 
The renormalization group equation which W[J] satisfies is, |llj . 

7m (g) 



9 or \ d 

V— + p(9)t^2 



dgi 



Jjj + [5(g)+2C lm (g)} — 



W[J] 



Thus from 
we find 

8(g) = 



6(g) 



2e + 2 lm (g) - 0(g) ^ 



(2.9) 

(2.10) 
(2.11) 
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The next stage, ^TQ, is to choose £(<?) to be the solution of the differential equation 



m d Tg ag) 



2im(g)((g) + 6(g) 



(2.12) 



(2.13) 



so that the running of the coupling coupling constant means that £(</) will also run according 
to its renormalization group equation, (|2.13|) . and which ensures that one can construct a con- 
sistent effective action which includes \ JA^ 2 as a composite operator. Given the nature of this 
differential equation its solution is singular at the origin and has the expansion, in powers of the 
coupling constant, 

oo 

£ c ng 2n • 



C(g) 



(2.14) 



n = -l 

With the explicit expressions we have quoted we can deduce that 
1 



g 2 C(g) 



(13C A - 8T F N f ) 
9N A 

+ (2685464C|T F iV / - 1391845Cj - 21M08C%C F T F N f - \mnmC A T F Nf 
+ 22ll84C A C F T F Nf + 584192C7 4 iV / 3 T| - 55296^7^/ 

i 4\ g 2 

- 65536T F N f ) 51Mjt 2 Na ( 35Ca _ i 6 r F JV / )(19C , A - 8T F N f ) 



'In the convention d = 4 — e was used. 
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+ ( (62228252520(7^ JV/Tf - 8324745975(71 - 42525100800(71 C^TV/Tp 

- 123805256256(7^ NfTp + 105262940160(71 C F iV/T| 
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+ 1 180127232(71 CpNfTp - 63001780224(71 <7 F iV/ 4 T F 

- 16782753792(71 NfTp - 475987968<7 A <7 F iV/ 4 T F 
+ 15308685312<7 A C F iV/T F + 3106406400<7aA/t| 
+ 70778880(7|iV/r| - 1478492160<7 F iV/T F 



- 2348810247V/ T F ) 99532^4^(35^ _ 16T pN f y(l9C A - 8T F N f f 
+ 0(g 6 ) (2.15) 

where the iV/ = case agrees exactly with the corrected result of the erratum of . 

Finally, for practical purposes it is more appropriate to transform the action, W[J], to a 
situation where it depends linearly on the source, J, which will allow one to apply the usual 
path integral machinery to construct the effective potential itself. In JTIj this was achieved by 
the Hubbard-Stratonovich transformation where one includes the redundant term 



1 = / Va exp ( - oio- + a 2 Ay + a 3 J ) (2.16) 



in the path integral where a is a scalar auxiliary field. The arbitrary coefficients, {a^}, are chosen 
appropriately to cancel off the \ J 2 and \ J A a ^ terms to leave a term proportional to o J as the 
only source dependence. These coefficients will depend on g, ((g), Z m and Z^. This results in 
a new Lagrangian, involving the field a, which includes the original QCD Lagrangian with its 
gauge fixing and a new set of interactions 

1 „n 1 /,-,/, ia\2 -aau n 1 aJM tAiM 



V 7 = - -G a uu G afiU (d^AlY - (fd^D.c 11 + i^ ll m i 

4 p 2a ^ 

jr—j— r — + — ^—aA a n A aix f^(A^A a A 2 (2.17) 

2gK{g)Z c 2g((g)Z c " 8((g)Z ( V * ) y > 



where 



e~ w[J] = J VA^VipVipVcVcVa exp J 



,r'.r I //' ^ 



(2.18) 



which is now linear in the source term. Here the quantity ((g) denotes the explicit expression 



1)2.15(1 and the extra terms again do not affect the renormalization of (|2.1|) which retains the 
usual renormalization constants. The sigma field corresponds to the composite operator from its 
equation of motion and is similar to the a field of the Gross-Neveu model where it is introduced 
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to remove the four fermion interaction. In (|2.17|l a new renormalization constant, Z^, appears 
which is defined to be 

z < = 1 + m <2J9) 

with £(g) also given by (|2.15|) . To the order we will require we note the value of the inverse, 
since that is what appears in (|2.17j) . is 
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6144^ 4 (35C A - lQT F N f ) 



We have now extended the LCO formalism of in the context of QCD and L°~ , (|2.17j) . 
has emerged as the effective action to compute the effective potential as a function of a. If 
a has a zero vacuum expectation value then one is in the usual perturbative QCD situation 
where the vacuum is known to be unstable, ^ [21 El 13- However, if a has a non-zero vacuum 
expectation value it generates a gluon mass in a non-perturbative fashion in a vacuum whose 
stability properties need to be checked. This requires the explicit form of the effective potential. 



3 Effective potential. 

With the previous formalism we are now in a position to construct the effective potential of 
the a field. It is derived by using standard techniques since the source for the operator is 
now implemented linearly in the generating functional. Beyond the tree approximation the 
calculation breaks into two pieces. Briefly the one loop correction is determined by summing a 
set of n-point diagrams with constant a external field. Hence, to one loop the potential is 
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9((g) 



(3.1) 



where the logarithm emerges from summing the appropriate diagrams based on ()2.17|) . However, 
the momentum integral is divergent and the infinity which cancels it arises from the counterterm 
available from Z^ in the tree term as indicated in our notation. Hence, expanding to the finite 
part one is left with the one loop effective potential 
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where we have set 
Ai = [13CU- 



8T F N f ] 



-1 



[3bC A -l%T F N f ]- 1 , A 3 = [19C A — STpNf]^ 1 (3.3) 
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9N A 



-a 



(3.4) 



(13CU - 8T F N f ) 

As a check on the result we note that (|3.2[) agrees with jllj in the Nf — > limit. To go 
to the next order we must piece in not only the higher loop Feynman diagrams but also the 
counterterms arising from the one loop integral of Q3.1|) . The latter is achieved by replacing 
&/{g({g)) in the logarithm in (|3.1|) by aZ m /(g((g)Z^). The resulting double and simple poles 
in e will cancel the same poles in the two loop corrections. These are determined by repeating 
the resummation of the one loop calculation with n-point a Green's function but now with one 
radiative correction included on each topology. However, it is well known that this procedure is 
equivalently reproduced by computing two loop vacuum bubble graphs constructed from a new 

This is given by shifting the a field to a new field, a, 



effective Lagrangian founded on 1)2.11 
which has a zero vacuum expectation value, and is defined by 



a 



(a) + a 



(3.5) 



and then dropping terms linear in a and the overall additive constant [44] . One, in principle 
does this for all fields in 1)2,17)1 . However, the gluon, ghost and quark fields all carry either a 
Lorentz or colour index which means that their vacuum expectation values are each zero and 
thus the shifted fields are equivalent to the original ones. The net result for 1)2.17)1 is that one 
has a massive gluon with a mass proportional to (a) and the interactions of this new action are 
effectively the same as before. Therefore, the two loop effective potential, in the Landau gauge, 
is deduced from 1)2.17)1 with a gluon whose propagator is 



1 



where 



(k 2 — m 2 ) 



m 



k»k v 
~k 2 ~ 



a 



(3.6) 



(3.7) 



In our case there are five graphs to compute. Two purely gluonic and three which involve 
either ghost, quark or a propagators. From 1)2.17)1 the propagator of the latter is momentum 
independent. 

The explicit calculation of each graph involves the evaluation of two loop vacuum Feynman 
diagrams with both massive and massless propagators. However, the exact expressions for such 
two loop topologies with three different masses is already known and given in, for example, (45] . 
We have written a Form programme to handle the substitution and determination to the finite 
part for each of the five diagrams. The result for the sum, after substituting for ((g) is 
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where S2 = (4\/3/3)CZ2(7r/3) and C^(^) is the Clausen function, ((n) is the Riemann zeta 
function and \x is the renormalization scale introduced to ensure the coupling constant remains 
dimensionless in ci-dimensions. It is related to the usual MS scale p, by p, 2 = Aire -1 ^ 2 where 7 is 
the Euler-Mascheroni constant. Extending the one loop calculation to the finite part of the next 
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order in the coupling constant and adding to (|3.8j) we find the finite two loop effective potential 
for non-zero Nf is 
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+ Cl 
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2048 



r 2, x2 , 3159 2 2 372015 
U^AiA 2 H — O pA\ A3 H — ; „„„ , A 



+ c5 



128 

+ ^C(3)C F A 1 A 2 
16 

_ 3 i^ C(3) c F A 3 2 
256 

/ 34749 



2048 

189295 
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8192 



V 256 



CfAiA 2 + 



1024 
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A^ 
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A 



*3 
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512 



1A2 



b4071 
1024 



A1A2 



lb384 

37179 w Nx , 37179 yl , 
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Q477 
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128 

12393 s , 2 37179 A/ ,, 2 
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409b 1024 
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C A + 
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409b 



Ca In 
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AiA 2 + 
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409b 
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2048 
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+ 
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1024 
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11 



7T J 



+ 0(g 6 



(3.9) 



which is the main result of this article. Again it agrees with the corrected potential of jllj when 
Nf = which is a strong check on our computations and renormalization though it is more 
complicated. 

With the explicit form we can make some general observations. One issue which always 
arises with potentials is that of boundedness. For a well defined physically interesting potential 
it must be bounded from below. However, the tree term is clearly dependent on the parameter 
Ai and therefore before any analysis can proceed we require 



A, 



[13C A -8T F N f ]- L > 



(3.10) 



as N A > 0. Interestingly this combination of group parameters is proportional to the one loop 
coefficient of the gluon wave function anomalous dimension in the Landau gauge. This would 
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suggest that it has to be positive for a bounded potential and clearly this cannot occur in the case 
of quantum electrodynamics, (QED). Moreover, there has been renormalization group studies of 
the relation of this particular coefficient of the gluon anomalous dimension in connection with 
confinement, |46l 147] . Thus, it is curious that a positivity condition emerges on the one loop 
Landau gauge anomalous dimension in the effective potential of a composite operator which has 
the dimensions of mass since the presence of a mass gap in QCD is indicative of confinement. The 
other main point about (|3.9j) is that we have computed the two loop corrections with massless 
quarks in the one two loop diagram where quarks can appear. As quarks are massive in the real 
world including a mass in the quark propagator in fact leads to divergent terms proportional to 
the quark mass which cannot be cancelled off by terms in the current effective potential. One 
might expect that a mechanism similar to that of the LCO technique could be used to introduce 
mass terms for the quarks. However, that approach runs into immediate difficulties since a four 
quark interaction, analogous to the four gluon term of ()2.17|) . would emerge which would clearly 
lead to a non-renormalizable effective action. 



4 Analysis. 



With the explicit two loop effective potential for arbitrary colour group we can now analyse 
various theories of interest. In particular we will focus on the SU(2) and SU(3) colour groups. 
As the latter is strictly the theory of quarks in the nucleons we will provide this analysis in 
detail. First, to appreciate how our analysis is carried out we consider the one loop potential for 
an arbitrary gauge group. From (|3.2j) the first task is to determine the stationary points defined 
by the solution of 

(4.1) 



which implies 



N A a' 



da 



9, /3. (ga'\ 3/13 207 \YI 

2 Al + U \W) + 128 + 1" 128 " ^° fX2 + ^ CfX " 

_ / 351 351 249 27 

+ Ca[ 3-CfAiA 2 + -t7-CfAiA 3 - — — A 2 + — A 3 

V 8 lb 128 b4 



+ c\ 



81. , 81, 
— AiA 2 + — A1A3 
lb 32 



(4.2) 



Setting the scale to be that given where the logarithm vanishes, ga' 
stationary points at 



a 







/i 2 , then there are 
(4.3) 



and 



where 



EL 

fjo 



= CAgl 36C A (35C A - lgZWVg) 

Vo ~ 16ir 2 (6545Cl-6008C A T F N f + 864C F T F N f + 1280Tl,Nf) 



(4.4) 



(4.5) 



This particular combination of the coupling constant on the left hand side has been chosen in 
order to compare with Clearly the trivial solution is the usual perturbative vacuum which 
has V(0) = whereas the latter induces a mass term for the gluon and is a global minimum 
provided (|3.1()j) is satisfied. With (|4.5|) we can obtain an Nf dependent estimate for the effective 
gluon mass, m e g, defined in ll a as 

m 2 eS = ga' . (4.6) 
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First, we need to convert our parameters to a reference mass scale, which is taken to be A-j^g 
and is introduced from the solution of the one loop QCD /3-function which determines the /j, 
dependence of the running coupling constant, g(fi), 



16vr 2 



2ft In 



AI- 
MS- 



where 



Hence, at one loop we obtain 



0o = y Ca 



A 



A {Nf] exD 
A MS 6XP 



(6545C^ - 6008C A T F N f + 864C F T F N f + 1280T|iVj 
2A(11C A ~ AT F Nf)(35C A - 16T F N f ) 



(4.7) 



(4. 



(4.9) 



which agrees with the one loop expression of jllj for the unitary colour groups when Nf = 0. 
It is worth noting that this formula is not valid in the QED case where Ca = 0, Cp = 1 and 
Tp = 1 because the original one loop potential is then unbounded from below. 

Repeating this analysis for the potential at two loops for an arbitrary colour group is ex- 
tremely cumbersome and it is more appropriate to concentrate on specific interesting cases. For 
instance, for SU(3) the two loop potential becomes 

36 



V(a)\ 



5(7(3) 




3537 



3537 



4(39-4iV / )(57-4iV / ) 
1851 



2(39 -4N f ) (105 -8N f ) 



8N f ) 



gV 2 

7T 2 
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1129005 



VL2 C(2) + ^C(2W 



+ 



128(39 -AN f ) (57 - AN f ) 

10161045 
512(39 - ANf) (57 - AN f ) 2 

1129005 

64(39 -4JV/) (105 -8AT/) 

10161045 



+ 



397305 



+ 



64(39 -ANf) (57- 4N f ) 
3575745 



256(39-4^(57-4^/) 
397305 



C(3) 
— C(3) 



32(39 -4iV/)(105- 8N f ] 
3575745 



128(39 -AN f ) (105 
3404293 

+ 



+ 



+ 



6144(57 -4^) 1024(57 -4^) 

397305 A/oX 7012085 
-C(3) + 



87V/) 2 ' 64(39 -AN f ) (105 - 8N f ) 

420267 1129005 

C(3) 



C(3) 
— C(3) 



1024(57 - 4iV / ) 2 
24525 



512(57 - 47V/) 
12222795 



6144(105 - 
132435 




32(105 - 8N f ) 



C(3) 



(4.10) 
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With (j4.1()j) we can determine the value of the coupling constant where the minimum occurs. 



N f 


„(i) 
Vo 


2 

Vo 





0.1925134 


0.1394790 


2 


0.2219195 


0.1832721 


3 


0.2398224 


0.2225696 



Table 1. One and two loop estimates of the coupling constant at the minimum of V(o~) for 

SU(3). 



These are recorded in Table 1 where we concentrate on the values Nf = 2 and 3 and the 
superscript here and elsewhere denotes the loop order. To consider more flavours one would 
be in the situation where the charm quark is present but for the case where it is treated as 
massless. As its mass is more comparable with the one loop effective gluon mass than the three 
light quarks and since we have neglected quark mass effects it is not clear how it would influence 
the approximations we have made. 

We can now determine the effect the two loop corrections have on the effective gluon mass, 
m e g . This requires the inclusion of the two loop correction to the solution of the /3-function for 
the coupling constant as a function of /x which gives 



16vr 2 



2/3 m 



A 2 — 

MS -I 



i -1 



01 



2/? 2 ln 



n -1 



A^ 

MS 



In 



2 In 



AI- 
MS- 



(4.11) 



where 



Pi = ~-c\ 



34 

y 



AC F T F Nf - ™C A T F N f . 



Converting the coupling constant at the minimum of the potential, go, when go' 



(4.12) 
jl 2 and 



N f 


m (1) /A W 


m (2) /A W 
m eff/ A MS 





2.03 


2.12 


2 


2.01 


1.99 


3 


2.00 


1.89 



Table 2. One and two loop estimates of the gluon effective mass for SU(3). 



substituting into ()4.6|) we obtain the one and two loop estimates given in Table 2. We have 

expressed our estimates in terms of the A^-dependent A^^ which for low Nf can be taken to be 
237 MeV giving a two loop estimate of around 450 MeV for Nf = 3. From the table, aside from 
the Yang-Mills result of jllj at two loops it is clear this effective gluon mass decreases slowly 
with the increase in the number of quark flavours. However, it is not clear whether this a real 
feature since quark mass effects, which should not be significant for light flavours, have been 
neglected. Moreover, it is worth stressing that m e Q is not the physical gluon mass in the new 
vacuum. The physical mass can only be deduced by computing the pole of the gluon propagator 
at two loops using the effective action (|2.17|) , This is clearly beyond the scope of the present 
article. 

We have repeated the above analysis for the case of SU{2) which follows the same path 
and the results are summarized in Tables 3 and 4. For iVy = 3 a second positive root emerged 

(2) 

for the value for which the potential was a minimum at two loops. This was Uq = 0.7559862 
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N f 




2 

Vo 





0.1925134 


0.1394790 


2 


0.2416107 


0.2239704 


3 


0.2758161 


0.4342486 



Table 3. One and two loop estimates of the coupling constant at the minimum of V(a) for 

SU{2). 



which we have neglected in the table as its value is significantly different from the one loop case. 
Though since the value given in Table 3 seems to be out of step with the pattern of Table 1 it 
is not clear how reliable that coupling is. The effective gluon mass values given in Table 4 are 
derived using the couplings of Table 3. Overall for most of the SU(2) and SU(3) cases we have 
analysed the two loop corrections to the effective gluon mass are no more than about 5% of the 
one loop result which is encouraging since it indicates a degree of stability of the approximation. 
However, the exception is the Nf = 3 SU (2) case which gives a gluon mass estimate which is 
20% different from the one loop estimate. As in this case a second positive root emerged for the 
coupling constant at the minimum, it may require a three loop effective potential to establish 
which if either is spurious. 



N f 


m (1) /A W 


m {2) /A W 
m cff/ A MS 





2.03 


2.12 


2 


1.99 


1.88 


3 


1.97 


1.58 



Table 4. One and two loop estimates of the gluon effective mass for SU(2). 

Finally, following the method of |31j . we have computed the value of the one and two loop 
effective potentials at the non-trivial minimum and presented the numerical values in Table 5 for 
SU(3) where they are all clearly negative. As the potentials are zero at the trivial minimum this 
justifies our remark that the non-trivial minimum corresponds to the stable vacuum. Clearly the 
inclusion of quarks increases the value of the minimum but the two loop correction is significantly 
larger. 



Nf 


mm/ V MS / 


v {2 ) I (a^\ 

mm/ V MS / 





- 0.32300 


- 0.76451 


2 


- 0.31145 


- 0.66966 


3 


- 0.30617 


- 0.61079 



Table 5. One and two loop estimates of the minimum of the effective potential for SU(3). 



5 Discussion. 

We conclude with various remarks. First, we have included quarks in the effective potential of 
the dimension two operator, \A a ^ in the Landau gauge and shown that the potential is bounded 
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from below for QCD but not for QED. Also, we have demonstrated that an effective gluon mass 
decreases in value slowly as Nf increases. However, this is for the unrealistic case of massless 
quarks, though for the values of Nt we have studied the quarks are light and their mass effects 
ought not to be significant. Clearly, for a more pragmatic study it would be better if quark masses 
could be included but this would require the removal of extra infinities which are quark mass 
dependent from the two loop effective potential. At present this will require an extension of the 
LCO formalism or another mechanism. Moreover, whilst our gluon effective mass estimates are 
essentially for a mass which from LF corresponds to a classical mass, a much better quantity to 
study would be the physical gluon mass itself which is defined as the pole of the gluon propagator 
derived from the effective action by computing the radiative corrections to the gluon two-point 
function. In the application of the LCO technique to the Gross-Neveu model and after improving 
convergence of the series, |32( I33| . it was shown that such a physical mass was accurate to a few 
percent with the known exact mass gap. It would be hoped that a similar calculation for QCD 
would give a reliable estimate for the gluon mass. Finally, given that we have now derived a 
two loop effective potential for a particular dimension two composite operator there is scope for 
examining the potentials of other related dimension two operators. For instance, ghost number 
breaking operators have been considered by various authors either in the maximal abelian gauge 
or in the Curci-Ferrari gauge which acts as a testbed for understanding gluon mass generating 
problems in the former gauge. Since we have simplified the way of computing the quantity 
without having to resort to the computation of three loop massive vacuum bubble graphs then 
the construction of such potentials for other operators or combination of operators can proceed 
in an efficient manner. 
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